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This paper considers the following generalization of the Weber plant location problem: the plant's 
output level is fixed, and its levels of input from its supply points, as well as its location, are among 
the decision variables. Hurter and Wendell (J. Reg. Sci., 1972) showed that this problem admits a 
kind of separability when the plant's production function lies in a certain class including the Cobb- 
Douglas forms. The present paper (a) determines the extent of that function-class, (b) carries out the 
explicit separation for the CES generalization of the Cobb-Douglas functions, and (c) discusses simple 
fixed-point-type iterative solution algorithms, similar to that well-known for the ordinary Weber problem, 
for several production functions (Cobb-Douglas, CES, and various two-stage technologies). Local 
convergence of these algorithms is established; computational experience will be reported in a separate 
Part II. 

Key words: CES; economics; Leontief; location theory; plant location; production functions; trans- 
portation; Weber problem, mathematical programming. 

1. Introduction 

The "ordinary" Weber plant-location problem, set in the real n-dimensional space /?", can be 
described as requiring the selection of xeR" to minimize the function 

m 

(f>w{x) = ^ t\\\x — Si\\qi. (1.1) 

i 

Here the decision variable x represents the location of a plant which requires m inputs for its 
operation; s, is the source of the ah input, t\ the associated unit transportation cost, and qi the 
level of the tth input. In (1.1), ||*|| denotes some appropriate norm on /?", which will be taken as 
the Euclidean norm throughout. The abbreviation 

pi= \\x-Si\\ (1.2) 

will prove convenient. 

This problem can be generalized by including the vector q= (q u . . ., q m ) of input levels among 
the decision variables. A plant output level q° is specified, as is the plant's production function 
f(q) and the unit prices pt of the inputs at their respective sources St. Now the problem is to choose 
xeR" and qeR m so as to achieve 

m 

minr min g ^ {hpi + Pi}q\ (1.3) 

i 

subject to f(q) = q°. (1.4) 
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A one-dimensional version of this problem was studied by Sakashita [1] l and extended to a net- 
work location problem by Wendell and Hurter [2]. The general case was formulated by Wendell 
and Hurter in a paper [3] which is the point of departure for the present work. 

The production function / will be said to lie in "class C(g )" if there is a positive constant 
K(q ) such that, for all positive q satisfying (1.4), 

m 
y Zq i dfldq i = K(q«). (1.5) 

It is noted in [3] that for this class of production functions, the problem has a kind of separability: 
it can be transformed into 

mm x /ji(x) (1.6) 

where /jl(x) is the Lagrange multiplier corresponding to constraint (1.4) for the "inner" minimization 
in (1.3). 

If / lies in class C(q°) for all q° > 0, we say it lies in "class C". As noted in [3], it is a con- 
sequence of Euler's theorem that class C contains all differentiate homogeneous functions, but 
it contains other functions as well. In section 2, we determine the extent of class C. 

Since class C contains the homogeneous production functions, it in particular includes the 
familiar Cobb-Douglas functions [4] 

f(q) = exp {a + J a < log*/,} (a* > 0) (1.7) 

as well as the multi-input CES ("constant elasticity of substitution") functions of Arrow, Chenery, 
Minhas and Solow [5], 

/(«)=(£ tar*)" 1 (fci>0,c>-l,c^0). (1.8) 

In [3] the separation (1.6) is carried out explicitly for the Cobb-Douglas case, but the corresponding 
problem for other homogeneous functions is noted to be "difficult". In section 3, we perform the 
explicit separation for the CES functions and for several functions representing two-stage tech- 
nologies. 

The ordinary Weber problem can be regarded as arising from a Leontief production function 
(in which specifying the output-level q° fixes the values qt), and is a convex programming problem. 
In contrast, the mathematical programming problems (1.6) arising from the cases studied here are 
nonconvex, so that their numerical solution is nontrivial. In section 4, simple iterative fixed-point- 
type solution algorithms are presented, patterned after one well known for the Weber problem. 
Local convergence is established for the low dimensions of practical interest, and the analyses 
necessary to handle the singularities s% are performed. A subsequent Part II will report our compu- 
tational experience, to date, with these algorithms. Further work should take up the case of time- 
varying prices p*, transport costs £*, and output levels q°. Another natural line of generalization 
would incorporate consideration of market location. 

2. Determination of Class C 

From (1.4) and (1.5), it is readily seen that C is the class of production functions /which sat- 
isfy a partial differential equation of the form 

m 

2<n9fldqi = F[f(q)] (2.1) 



1 Figures in brackets indicate the literature references at the end of this paper. 
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for all positive q. Since a production function / has positive first-order derivatives, F is positive- 
valued; we also assume it continuous. 

THEOREM: Class C consists of the production functions of the form f(q) = M[h(q)], where h is a pro- 
duction function homogeneous of degree 1, and M is an increasing differentiate function. 

Before proving this theorem, we note three consequences of it. First, it identifies class C with 
the class of "nomothetic" production functions introduced by Shephard (see p. 30 of [6J), apart 
from questions of smoothness and other properties in a theoretical definition of "production func- 
tion". Second, it "explains" the examples of nonhomogeneous members of C given in [3], which 
in fact are the logarithms of Cobb-Douglas functions. Third, it implies that analyses of our gener- 
alized Weber problem can be confined to production functions which are homogeneous of degree 1, 
since with /= M[h\ as above, the constraint (1.4) can be replaced by the equivalent h(q) = M~ x (q {) ). 

For the proof, first assume f=M[h] as in the theorem's statement. By Euler's theorem on 

homogeneous functions, m 

^qidh/dqi = h(q), 
i 

and so by the chain rule of differentiation 

in 

^qidfldqi = M'\h{q)]h{q), 

1 

an instance of (2.1) with F(v) =M f [M" 1 (v)]M~ l (v). 

Conversely, suppose/eC satisfies (2.1). Define a function M~ l (u) by 

M-*(u) = exvl [" [l/F(v)]di!\ (2.2) 

for some u {) > 0. Since F is positive-valued and continuous, M~ x is increasing and differentiable, 
and thus has an increasing differentiable inverse function M. Define h — M' l [f]; then by (2.1) 
and (2.2) 

in m 

£ qidhldq,= (M -')'{/] ^qrffldq, 

1 1 

= {M->[f]/F[f]}F[f] = h, 

so that the converse of Euler's theorem implies that h is homogeneous of degree 1. Since /= M[h] , 
the proof is complete. 

Before leaving this topic, we note that class C also arises in a multi-output generalization of 
the problem under discussion. Namely, suppose the m inputs are used jointly to produce several 
outputs in accordance with a vector production function having one component fj(q) per output. 
Suppose furthermore that the level of each output is prescribed, and that p,j(x) denotes the 
Lagrange multiplier corresponding to the jth of these constraints in the inner minimization of 
(1.3). The first-order optimality conditions for that minimization are 

tiPi+Pi=^LW(xWjldqi (i=l, 2, . . . ,m), 

3 

so that with the inner minimization accomplished for each x (and the dfj/dqi evaluated at its solu- 
tion) the objective function (1.2) becomes 

^ {tipi + Pi)qi=^ Vj(x) 2 QidfjIdQi- 

i J I 
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If now each fj is in class C, with corresponding Fj in (2.1), then this minimand is equal to 

4>(x) = ^ H (x)Fj(ff) (2.3) 

J 

where ff is the prescribed level of the jth output. Result (2.3) is a multi-output extension of the 
separability expressed by (1.6). This (multiple output) line of generalization will not be pursued 
further in the present paper. 

3. Some Cases of Explicit Separation 

To take full advantage of the separability expressed in (1.6), it is necessary to find an explicit 
expression for (Ji{x) , so that the resultant "pure location" problem is in explicit form. This will be 
possible, in particular, if/ satisfies a suitable set of identities 

q i = F i [f(q),df/dq i ]. (3.1) 

To see why this is so, recall the first-order optimality conditions 

tipi+ pi= fJi(x)dfldqi 

for the inner minimization in (1.3). The solution q(x) of that minimization will therefore satisfy, 
by (3.1), 

qi(x) = Gi[(tiPi + pi)lfi(x)] 9 

where Gi(')=Fi(q° 9 -). It follows that 

q° = f[q(x)]=f{Gi[(t l pi+p l )ltJL(x)],. . . 9 Gn[(tmpm+Pm)lli(x)]}- (3.2) 

Typically (hence the adjective "suitable" above (3.1)) this equation can be solved to obtain the de- 
sired explicit form for p(x). 

For the Cobb-Douglas case (1.7), one has in (3.1) 

Fi(u, v) = aiu/v; 
thus (3.2) yields 

m 

log q°=a+ 2 <H log[aiq°p(x)l(tipi+Pi)] 

i 

m m m 

= [a+]T a,- log (a,<7°)]+ [logfi(x)] £ ai-^at log (tipi+pi)' 

i i i 

It follows that 

/ m \ - 1 m 

log p(x)= \Y at J ^ ai lo S (*iPi + Pi) + const. 

^ i ' i 

or, with the abbreviation 



7Ti = pj/ti, 

)-l m 
2 ai log (pt+7Ti) + const. 
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(3.3) 



Thus the pure location problem (1.6) is equivalent in the Cobb-Douglas case to minimizing 

<J)Cd(x) = ^ di log (p/+7T/), (3.4) 

1 

a result derived in [3|. 

For the CES production function (1.8), one has in (3.1) 

Fi(u, v) = u(bilv) 1 f< c+1 \ 

With d=c/(c -hi), (3.2) yields 

III 

(? )- f =^i[? ,, (^(^/{^+ft!)" lr+l, ]- r 

1 

m 

= (q»)-r[fx(x)]-«^b)-«{t i p i + p l }" 

= «)-'[ M (^)]-"£ Ci {p i +77,}" 



where ir,- is as above and 



Ci=b\- d tf. (3.5) 



If c>0, so that 0<rf<l, then minimizing /jl(x) is equivalent to minimizing [/x(x)] r/ , and it 
follows that the pure location problem for the CES case can be expressed as demanding the mini- 
mization of 

m 
<£r/;.v(*) = £ C,{p, + 77-,}" (0<rf<l). (3.6) 

1 

If d < (i.e.,— 1 < c < 0), the pure location problem involves maximizing the form (3.6), or 
equivalently minimizing its negative. Our subsequent discussion of the CES situation is readily 
adapted to this subcase, but will for simplicity be confined to the subcase < d< 1; the reader is 
warned that the later discussion does not as it stands refer to the case c < (i.e., d < 0), though the 
revision is simple. 

The three functions c/>» , <{)cd, 4>ces all have the form 

m 
cM*) = £ <fc(pi) (3-7) 

1 

where the functions <pi are defined and nonnegative on (0, oo) 9 positive-valued and twice differen- 
tiable on (0, o°), and satisfy 

4>;>0 on (0,oo). (3. 8 ) 

But while for (\>w> which has <f)i(u)= (ttqi)u, each summand in (3.7) is a convex function of x, 
neither $r/; nor (J)ces is convex, so that the CD and CES cases give rise to nonconvex programming 
problems. This nonconvexity is most easily seen in the one-dimensional case; whereas (j)w is 
linear on each (open) interval between successive points 5/, both (/v/; and (f)(Ks satisfy <f>" < (the 
antithesis of convexity) on those intervals. 

The absence of convexity suggests the possibility of multiple local minima, and these can in 
fact occur. They may occur at a point s, (in the one-dimensional case, local minima occur only 
at points si), which however would not be routinely identified as a critical point since Si is a singular 
point of 
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gradp,= (x — St) I pi 

and thus of grad </>. We therefore proceed to develop a special test for the existence of a local 
minimum at an s;, saysi. It will be assumed that 

Si = Si for 1 ^ i ^ r, Si^Si for i > r. 

Let be a nonnegative scalar variable, and weR n be a variable "direction vector", i.e., || w || 
= 1. Set 

g(0,w)=<l>(s 1 + ew). (3.9) 

Then a necessary condition for a local minimum at s i is that 

M w dg(0+ 9 w)ld0&0 (3.10) 

holds. For x = s i + Ow, one has 

Pi =6 for i ^ r, p,-= ||0w;— (s,- — Si)|| for i > r, 

so that 

/• m 

g(6,u>) = 2 *<(») + £ 4>i(\\0w- (si-s,)||), 

1 r+1 

a^/ae = J *»' < e) + 2 [*i'(p«)/p«] [»- (»a-*i)] (3.11) 

1 r+ 1 

where (w,s; — Si) denotes the scalar product. With the notations 

m 

b= £ fof(||*<-*i||)/||»i-ai||](*-4i). (3-12) 

r + l 

A = j^<j ) ' i (0), (3.13) 

1 

it follows from (3.11) that 

dg(0+ 9 w)ld6 = A - (w, b). (3.14) 

The Cauchy-Schwarz inequality implies 

(w,b)^\\b\\, 

with equality for all w if 6 = 0, and for w=6/||6|| and its negative if b ^ 0. It follows from (3.14) that 

min w dg(0+, w)dO = A-\\b\\ , 
and so 

A&\\ b || (3.15) 

is a necessary condition for a local minimum at si. 
Conversely, suppose 

<r = A-\\b\\>0. (3.16) 

Choose any positive 8, < min;> r | |s ,-s; 1 1. Then the right-hand side of (3.11), call itg*(0, w) , 
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is continuous on the compact domain [0, 81] X {w: || w\\= 1}, hence uniformly continuous. In par- 
ticular there is a 6 > 0, with 8 ^ 8 U such that if ^ 6 ^ 8 then for all w 

g*(6,w) >g*(0,w) -o-. 

Since the previous analysis implies #*(0, w) ^ cr, we have demonstrated the existence of a 8 > 
such that 

dg(0, w )ldd=g*(0, w)>0 for ^ ^ 8 and all w. 

Thus in the 8-neighborhood of Si, </> is uniquely minimized at the point 5i, so that (3.16) is a suffi- 
cient condition for a strict local minimum at si. It is a generalization of its specialization (given by 
Kuhn and Kuenne [7] ) for the ordinary Weber problem. 

Next we consider some cases in which the "inputs" transported to the plant from the source- 
points 5, are best interpreted as "factors of production", not for the process yielding the plant's 
final output, but rather for intermediate on-site processes producing these "final factors." Note 
that the production functions for these intermediate processes, as well as that for the final process, 
must now be specified. The levels of the final factors are (intermediate) variables of the problem; 
these levels will be denoted 

Q=(Q„ <?., . . .,£>») 

and the production function for the final process will be denoted f(Q). 

A variety of interesting problems can be posed in this context; we will briefly take up just a 
few of them. For notation, it will be convenient to partition the input-indexing set {1,2, . . . , m} 
into subsets {/(i>):i>=1, 2, . . ., A/}, where iel(v) signifies that the ith input goes into making 
the ^th final factor. Assuming disjointness of these sets I(v) is not really a restriction on the tech- 
nology—so long as capacity constraints at the sources are omitted — since otherwise-identical 
inputs can be artificially distinguished according to the final factor in which they will be embodied. 

Suppose first that each intermediate process follows a simple Leontief production law; that 
is, there are positive constants Ki such that 



qi = KiQ v foralli€/(i/). (3.17) 



Then the problem can be written 
subject to 



minj. min Q ^ ^ {tip\ + p,-} K\ 

iel(v) 



0, (3.18) 



f(Q) = q°. (3.19) 



Reduction to a pure location problem follows the same pattern as before; if/ is a Cobb-Douglas 
function with parameters a v , or a CES function with parameters b v and c, the result is an objective 
function 



<\>'co (*) = 1 a " lo g £ MP* + in] ) , (3.20) 

UHp) 



nth Li = Kiti, TTi= pilti, or 



Vces (*) = I V-" ( 1 Kdhpi + Pi} )" 

" \ hld>) ' 
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with d= c/(c + 1), 7T; as above, and Li = b l J c Kitt. 

Next suppose that the final process is of Cobb-Douglas type, with parameters a v . If each 
intermediate process is also Cobb-Douglas, then the qt's are related to the plant output by a Cobb- 
Douglas function, so that the material leading to (3.4) applies. Let us suppose, instead, that each 
intermediate process is of CES type; assume the ^th final factor has aCES production function with 
parameters {bi'.iel(p)} and c„. The result is a (composite) production function, for the plant, 
of Uzawa-CES type [8], 



ift(q) = exp la— 2 {aJc v ) log ( ^ Mr^jf • 

It is readily verified that 

dilfldqi=a v biqi( c v +1) \l* / ]£ b } q- c » (iel(v)) 



(3.22) 



(3.23) 



je/(i>) 



Although identities of type (3.1) are lacking, the general approach can still be carried out. Let 

ai(x) = tipi + pu (3.24) 

then (3.23) and the first-order optimality conditions below (3.1) yield 



a t = navbiqr^+Vif, / 2 hqfv (iel(v)), 

or equivalently, with d v = c v \(c v + 1) , 

q -c v = {^)-d v { ai la v bi)-d> ( 2 bjqreuV (iel(v)). 

Multiply both sides by bi and sum over iel(v); the result is 



iel(i') \jel(v) ' iel(i') 



or equivalently 



2 biqr c v= {^a v )- c v ^ biiajbi)^ 



Ul{v) 



iel{v) 



<V + 1 



i.e., 



log (^ b iq ^) =-cJog ft- cAog(iPa,,) + ( Cl ,+ l) log £ 6/(^/^)^|, 
leading via (3.22) to 

iog<f =a+( 2 a log ^" h 2 a ^ log ^ ()a ^~2( a ^) log 2 W^M*']- 

Thus, with the abbreviations 

$ v = a v \d v , 7Ti = pilti,Li=b l .- d vtfv, 
the objective function for the pure location problem takes the form 

4>£f(*) = 2> log 2 £<{pi+^} d " 



(3.25) 
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If the/„'s are as above, but /is a CES function with parameters b v and c, then the composite 
production function is 



Uq)= 2 &, XK 



i€l(v) 



elc u -l/c 



Manipulations like those above lead to an objective function 



''/"„ 



(3.26) 



(3.27) 



with d=c/(c + l), rfv=c„/(c„+l), 7Ti=pilti and Li = b*v fc b}- d »tf'>. 

Now assume /is as above, but the/Vs are Cobb-Douglas functions with parameters a„ and 
{a,-:ie/(i^)}. Then the composite production function is 



*/>(<?) = ^b v exp (a^ + ^ at log g< 



I €/(!/) 



-c L-l/c 



(3.28) 



For the explicit separation to be tractable, it appears necessary to require each/, to be homogeneous 
of the same degree, i.e., 



£ a i = ^ (all v). 



(3.29) 



With the notations 



an objective function 



i€l{v) 



cr = c/(l + cA), A v =a u + V a/loga,, 

i€/(v) 

7Ti = pi/bi, a\ = (Tai 

a\= ^ a! log ti+ (1 — o-A) log6„ — av4„, 



0c2sW = 2 exp | a J + 2 a i log (pi + ir,) fi 



(3.30) 



a sum of Cobb-Douglas functions, is obtained for the pure location problem. 

The reasonableness of the restriction (3.29) is supported by the following observation, which 
applies to the situations (3.25), (3.27), (3.30) above. Suppose /eC, with associated function F in 
(3.1), and that the functions f v are homogeneous of respective degrees A„. For the two-stage tech- 
nology to admit the kind of analysis given in this paper, the composite production function i// must 
lie in C. Now 

2*#/«* = 2d//^ 2 qidfjdqt 
= ^ A v Q v dfldQ v , 

v 

and only if all A„ have a common value A can we continue to the 

= AFM 

which shows that \\)eC. 

We return now to the matter of testing for a local minimum at a point su Since the objective 
functions (3.20), (3.21), (3.25), (3.27) and (3.30) are not of the form (3.7), the test (3.16) does not 
apply. Instead, these objective functions have the more general form 
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000 =^<t>Api>) (3.31) 

V 

where p v is a vector with components {pi\iel(v)}. For the point 5 being tested, define the index- 
sets 

E(v) = {iel(v):si = s}, U{v) = I(v) -E(v), 

and introduce a variable vector u v with components {ui:iel{v)} as general argument of <f> p . As in 
the analysis leading to (3.16), let 

g(0, w) = (f>(s + Ow) 
with 6 a nonnegative scalar variable and weR" a direction vector. Then 

i> iel{v) 

[6-(w,Si-s)]l\\dw-(si-s)\\ 
= 2 2 [«*»/*»«] [{||«w-(*,-«)||W)] 

+ E 2 [^/3ui]({||0«;-(^- s)||}>#w] 

k ie£7(i>) 

[0- (»,«|-*)]/||ft«7- (Si-S)||. 

It follows that 

d*(0+,«;)/afl = 2 X [^"iKflk-'IIW 

~2 E [ a ^/^]({|l 5 J _5 llbc/(*'))(^ 5/-5)/||5/-s||. 

Arguing as in the derivation of (3.16), we obtain the criterion 



where now 



A>\\b\\ (3.32) 

/l = IX [d<t>JSui]({\\si-s\\} j(IM ), (3.33) 

6= 2 2 (*-*)[^/a«N]({||*j-*||}i.«r))/||*»-«||. (3.34) 

Note that in (3.33) the arguments \\sj — s\\=0 for jeE(v). 

One might also consider a two-stage process with the final stage of Leontief type. But then 
fixing q {) fixes all Q v , so that the problem is equivalent to one of the single-stage multi-output type 
described at the end of section 2. 

The final situation to be considered is that the various inputs iel(v) are of the v\h final factor 
itself (without further processing), but are distinguished merely by being from different sources. 
That is, the v\h intermediate process has as "production function" the additive 

Q v = 2 *• < 3 - 35 ) 

iel(v) 

Here formal use of the preceding approach would lead to nonsense. The reason lies in the reliance 
of that approach upon the optimality condition below (3.1) to characterize the inner minimum in 
(1.3). In fact, that condition is guaranteed only for those qi which are strictly positive at the mini- 
mum, a condition which indeed is satisfied in all the previous cases treated, but is violated here 
since each final factor would be purchased only from the least expensive of its sources. 
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With ai(x) defined as in (3.24), let 

a*(*) = min {a t (x) :iel(v)}; (3.36) 

then (1.3) for the case (3.35) can be written 

min min ^T a*(x)Q t ,. 

x Q v 

If for example /is Cobb-Douglas with parameters a v , the resultant pure location problem has in 
analogy with (3.4) the objective function 

*?*(*) = 2 fl * lo « [«*(*)]• (3.37) 

V 

while if /is CES with parameters 6^ and c, the result is 

CW = S^KW] d (3.38) 

where d=c/(c+l), analogous with (3.6). Under the plausible assumption (in the present context) 
that {ti \iel(v)} has a single member £„, (3.37) and (3.38) can be replaced by 

,*„(*) =2 a " log ^ min {p* +7r *}]' ( 3 - 39 ) 

**«(*) -2«» [ min {P«+Wi}]" (3.40) 

where c v =b\~H { l and 7T\ = p\\t\. Note that under the further assumption that {pniel(v)} has a 
single member p„, 

min {p,-+ 7T,-}= [min p,-] + 7r t , 

t€/(v) 16/(10 

with 7T v = Pvlt v . 

4. Iterative Solution Methods 

The pure location problems obtained in section 3, by working out several cases of "explicit 
separation", require the minimization of fairly complex nonconvex functions (f)(x). Since such 
problems are computationally nontrivial, it seems useful to present a class of iterative solution 
methods which are simple in concept and simple to program. These algorithms, which are based 
on characterizing an optimal location as a fixed point of an associated transformation of/?", are 
presented in the present section, while computational experience with them will be reported in 
Part II. 

As noted in (3.7) and (3.8), several of these problems have an objective function of the form 

m 

</>(*) = ]£ <MP/) (4.1) 

i 

where the functions 4>i(u) are twice differentiate for nonnegative arguments, and satisfy 

<$>\ >0 (£ = 1,2, . . . ,ro). (4.2) 
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At any point x not an s i , one has 



grad <f)= ]£ (/>; (pi) (x - 5t)/pi- 

1 

If x is to be a local minimum, then grad 0=0 must hold, or equivalently 

m 1 / r 

^Si&ipdlpi / 5>,' (p«)/p 



(4.3) 



(4.4) 



a formula which displays x as a fixed point of the function on the right-hand side, and incidentally 
as lying in the convex hull of the points 5/. This formula suggests the iterative scheme 



Ak + \) — 



2 5^; (pf)ipf 



2 *; (pf)ipf 



(4.5) 



Because of the presence of (in general, uncancelled) denominators p\ k \ this form can be unsuitable 
for numerical work when x {k) is near some Su say Sj, and should be replaced by the algebraically 
equivalent form obtained by multiplying numerator and denominator through by pj. A) . (This alterna- 
tive form also shows that each Sj is a fixed point of the transformation given by (4.5).) Of course, 
the test (3.16) for a local minimum at Sj should be applied in such cases. 

For the ordinary Weber problem, with objective function <f> w given by (1.1), the algorithm reads 



r (fc+i) = 



X SitiQilPi 



(k) 



2 UqM 



(k) 



(4.6) 



This iterative scheme, which has been repeatedly rediscovered (e.g. [7], [9-111), goes back at 
least as far as Weiszfeld [12], who also gave a convergence proof; the rapidity of that convergence 
has been confirmed in a number of instances, e.g. [13]. 

For the Cobb-Douglas case, with objective function $ CD given by (3.4), the algorithm reads 



Ak+D — 



£w(pi*>+iri)pi* 



^aiKpW + ndti 



(k) 



while for the CES case, with objective function </> given by (3.6), it reads 



Ak + D = 



JsiCi/Cp^ + TTi) 1 "^ 



^cJipP + iriV-W 



(4.7) 



(4.8) 



This scheme (4.8) was considered by Cooper [14] for the case of all 7T ? = 0, a limiting case of the 
situations of interest here. Note that if some 7Tj = 0, and if the algorithm leads to an x (A) near 5j, 
then the numerator and denominator need to be multiplied by [p^] 2 in (4.7), and [p { j ) ] 2 ~ d in 
(4.8), not just p^ A) . Note also that if some 7Tj=0, then for both <f> and <j) the test for a local 
minimum at Sj yields A = o° in (3.13) and thus an affirmative result for the test; for , with all 
ttj = 0, the fact that each Sj yields a local minimum was observed by Cooper [15]. 

By the local convergence property (LCP) for the pure location problem, we shall mean that each 
strict local minimum x of (/), other than the points s,, has a neighborhood N(x) such that if the itera- 
tive process enters N(x) at some stage, then it subsequently converges to x (in fact, in an-least- 
geometric fashion). In a paper [16] dealing with the general scheme (4.5), Katz (op. cit., Theorem 4) 
shows that the LCP holds if, in addition to (4.2), the functions (/), satisfy 



4>;»^(3-/i)(/>;»/u. 
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(4.9) 



For both <b , and <£>,,,, one has </>" < 0, so that for the low dimensions (n ^ 3) of greatest practical 

(I) ( h & ' 

interest, (4.9) is satisfied and hence the LCP is assured. 

The objective functions (3.20), (3.21), (3.25), (3.27) and (3.30) have the more general form (3.31), 

<M*) = 2 4>»(Pp), (4.10) 

with p v the vector with components {pi\iel(v)}. The twice-differentiable positive-valued functions 
<\) v (u v ) , where u v is a vector of nonnegative variables {u,:iel {v)} , satisfy in all these cases the 
analog 

d<f) v ldui>0 (aHid(v)) (4.11) 

of (4.2). The analysis by Katz [16 1 can be mimicked to obtain a generalization of (4.9) which, together 
with (4.11), is sufficient to assure the LCP for the generalization (given later, below) of (4.5). 

The details of this imitative analysis are straightforward by reference to [16|, and therefore will 
not be repeated here. The result is that the condition 

'* I iel(v) i,jel(i>) ' 

(4.12) 

together with (4.11), sufficies for local convergence at x. It follows that the conditions 

(n-3) ^ (l/p«)[a^/3tt|](p„)+ ^ [(x-s u X-sj)lp i pj][a^ v ldutau J ']{p v ) ^0 

iel(r) i,jel(i>) (4.13) 

for all v, together with (4.11), are sufficient. In particular, if for each v the local minimum x lies 
outside the convex hull of the points {s;:ie/(*>)}, so that in (4.13) each scalar product (x — 5t, 
x — Sj) > 0, and if each 0„ has all d 2 (f)v/diijdiij ^ 0, and if n ^ 3, then local convergence holds at 
x. For a more useful condition, one can employ the consequence 

(x — Si, x — Sj)lpipj*2 (-1) (4.14) 

of the Cauchy-Schwartz inequality. If each $ r satisfies 

d^JdUidUj^O (ij€l(v); i^j), (4.15) 

then it follows from (4.13) and (4.14) that 

^ {(n-3) {llui)d<t> v ldui + d 2 <£„/dttf} -^ {d 2 <f>vlduiduj:ij€l(v);i^j}^0, (4 16) 

together with (4.11), is sufficient to assure that LCP. Note that (4.16) is a generalization of (4.9). 

Consider first the objective functions (j)^ D and (f)^ ES of (3.20) and (3.21). For each of them, 
(f> v (u v ) has the form 

4> v (u v ) = g v ( 2 LiWi + m}\ (4.17) 

so that (4.16) takes the form 

X {(n-3) (LJmW + L'gl'} -^{L i L j g' l ,':i,jd(p),i^j}^0 (4.18) 

iel(r) 

with gi and g" evaluated at the g„-argument of (4.17). For (3.21), with gu(v)=v d , this condition is 
^ \(n-3) (Ltlui) ^ Lj(uj+7rj)-(l-d)q\ + (1 - d) £ {LiLj:iJeI(v); i #/} ^ 0, 
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or equivalently 
£ Lf{(n-3) (Hut) (m + irO-d-rf)} 

+ 2 {Wi[(»-3){(u j + ff j )/B,+ (ai + ir*)/^} +.2(l-rf)]:i,jtf(y),i</}^0. 

Since rf < 1, assuming n ^ 3 assures that the first sum is < 0. The generic summand of the second 
sum, divided by LiLj, is (if n ^ 3) 

^ (/i-3)[a i /i4«+tti/ttj]+2(l— rf); 

applying to the first term the inequality z+ II z ^ 2 forz > 0, proves for n ^ 3 that the last displayed 
expression is 

^2(n-3)+2(l-d)=2(n-2-d), 

which is negative for n ^ 2. Thus, for <f)^ ES , LCP holds for the planar and one-dimensional cases. 
The same argument, with d=0 in the later steps, yields the same conclusion for <f)^ D . 

For the objective functions 4> c ( ^ and 4> C C E E S S of (3.25) and (3.27), we have the generalization 

(f} v (u v )=gJ 2 Liim + in}**) (4.19) 

of (4.17). Thus (4.16) takes the form 

2 L i (u i ^7T i ) (i ^[{(n-3)(u i ^7T i )lu i -(l-d l/ )gl^L i (u i +7T i ) (i ,gldJ 

i€l(v) 

-d p ^ {LiLj (m + 77/ ) d v~ x (uj+TTj ) d v' x g'l : i, jel(v) ;i^j}^0. 

For (3.27), vnthg v (v) = v 8 v where b v = dld v , this yields 
^ L i (u i ^7T i ) d ^[{(n-3)(u i 7r i )lu i -(l-d l> )} X L^+irj)^ 

iel(v) jtl(v) 

+ L i (u i +Tr i )<> l ,(d-d v )]-(d-d„)"2{L i Lj(u i +TT i )< i „- l (u j +iTj)\- , :i,jeI(v);i^j}^0 
or equivalently 

£ Lf (m + TTi) 2 ^- 2 [(n - 3) ( Ui + TTi)lm - (1 - d v ) + (d- d v )] 
hiM 

+ J j {L i L j (u i + 7T i )<>^ (uj+ir,)**-* [{(i»-3) (ui + ir,)/u«-(l -<*,)}(«, + «■,,)* 
+ {(n-3) (» J + ir,)/B j .-(l-rf,)}(ui + ir 1 ) I -2(rf-rf,) (u« + in) (bj+it,)]} 

:i,jel(v);i<j}^0. 

Since a*<l, the first sum is negative for n =£ 3. As for the second sum, the factor [ ] in its 

generic summand is forn =£ 3 

«-[{(l-^)-(n-3)}{(u i +77 i ) 2 + (Uj+Tr j y i } + 2(d-d„)(u i + TT i )(u J +1Tj)] 

= -{(l-d,)- (n-3)}<?(«i+ ir f ,itj+irj), 
where the coefficient (1 — d v ) — (re — 3) is positive for n =£ 3, and the quadratic form Q is given by 

Q(v,w) = v 2 + w 2 + 2kvw, k= (d-d„)l{(l-d„) - (re-3)}. 

Q is positive definite for k 2 < 1. Since k < 1 follows from the fact that d < 1, it suffices to have 
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k> (—1). For n = 3 this is true if d + 1 > 2d„, while for n ^ 2 it follows without additional restric- 
tion. For (3.25) the analysis is similar, corresponding to d = 0. 

For the objective function <f>cEs °f (3.30), condition (4.16) leads to 

2 a\ (in+ 7n)~ 2 {{n-3) (ui + 7Ti)/ui — l} 

icl(v) 

which holds for n ^ 3. 

The generalization of (4.5) to the situation (4.7), for which the preceding convergence analyses 
employing (4.16) were given, is obtained using the generalization 

grad = 22 \- d< f>*l du i] (P") (* - s i)lpi 



of (4.3). The iterative scheme is 



r (A + D = 



2 2 Si[ d <t>vldUi](pv)IPi 
v iel(i') 



2 2 [^/^,](p,)/ P , 



(4.20) 



For x (k) near some s,-, precautions like those noted under (4.3) are required. 

It remains to consider objective functions like the c/>* /; and <f)* FS of (3.37) and (3.38), whose 
general form is 

<f>*(x) = 2 <£i'[«*U)L a*(*) = min {<*,-(*) :i'e/(^)} 

with a, defined by (3.24). Let us call x exceptional if a tie occurs in the definition of some a*(x); 
for any nonexceptional point x, let i(v, ^)e/(^) be the unique index for which the minimum occurs. 
For nonexceptional points, the previous analyses can be carried over by replacing a*(x) with 
<Xi( l ,,x)(x)> The reason is that these analyses — testing for a local minimum at a points/, or for local 
convergence to a local minimum (not an Si) of an iterative scheme based on a zero gradient (itself 
a local construct) — deal only with local behavior of $*, and each nonexceptional point x has a 
neighborhood consisting entirely of nonexceptional points y for which i(v, y) = i{v, x) for all v. 
But the "radius of convergence" around a local minimum is reduced by the need to avoid contact 
with the set of exceptional points; a local minimum lying near this set may therefore be "hard to 
get at" for the algorithm. If the algorithm generates an x {k) which is an exceptional point, it is 
natural to proceed by breaking the tie arbitrarily, and the effect of this seems difficult to predict. 
(An alternative is to employ a more complex logic involving "branching" when an exceptional 
jc (a) is encountered.) The ability to detect a minimizing point which is exceptional is a priori dubi- 
ous. These problems are explored on an empirical basis in some of the computational experiments 
to be reported in Part II. 
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